A new model of binary opinion dynamcis: coarsening and effect of disorder 
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We propose a model of binary opinion in which the opinion of the individuals change according to 
the state of their neighbouring domains. If the neighbouring domains have opposite opinions, then 
the opinion of the domain with the larger size is followed. Starting from a random configuration, the 
system evolves to a homogeneous state. The dynamical evolution show novel scaling behaviour with 
the persistence exponent 6 ~ 0.235 and dynamic exponent z ~ 1.02 ± 0.02. Introducing disorder 
through a parameter called rigidity coefficient p (probability that people are completely rigid and 
never change their opinion), the transition to a heterogeneous society at p = 0^ is obtained. Close 
to p = 0, the equilibrium values of the dynamic variables show power law scaling behaviour with p. 
We also discuss the effect of having both quenched and annealed disorder in the system. 



PACS numbers: 87.23.Ge, 68.43Jk, 87.15.Zg, 89.75.Da 



A number of models which simulate the formulation 
of opinion in a social system have been proposed in the 
physics literature recently [H] . Many of these show a close 
connection to familiar models of statistical physics, e.g., 
the Ising and the Potts models. Different kinds of phase 
transitions have also been observed in these models by 
introducing suitable parameters. One such phase transi- 
tion can be from a homogeneous society where everyone 
has the same opinion to a heterogeneous one with mixed 
opinions 

In a model of opinion dynamics, the key feature is the 
interaction of the individuals. Usually, in all the mod- 
els, it is assumed that an indiviual is influenced by its 
nearest neighbours. In this paper we propose a one di- 
mensional model of binary opinion in which the dynamics 
is dependent on the size of the neighbouring domains as 
well. Here an individual changes his/her opinion in two 
situations: first when the two neighbouring domains have 
opposite polarity, and in this case the individual simply 
follows the opinion of the neighbouring domain with the 
larger size. This case may arise only when the individ- 
ual is at the boundary of the two domains. An individual 
also changes his/her opinion when both the neighbouring 
domains have an opinion which opposes his/her original 
opinion, i.e., the individual is sandwiched between two 
domains of same polarity. It may be noted that for the 
second case, the size of the neighbouring domains is ir- 
relevant. 

This model, henceforth referred to as Model I, can 
be represented by a system of Ising spins where the up 
and down states correspond to the two possible opinions. 
The two rules followed in the dynamical evolution in the 
equivalent spin model are shown schematically in Fig. [1] 
as case I and II. In the first case the spins representing 
individuals at the boundary between two domains will 
choose the opinion of the left side domain (as it is larger 
in size). For the second case the down spin flanked by 
two neighbouring up spins will flip. 

The main idea in Model I is that the size of a domain 
represents a quantity analogous to 'social pressure' which 
is expected to be proportional to the number of people 
supporting a cause. An individual, sitting at the domain 
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FIG. 1: (Color online) Dynamical rules for Model I: in both 
cases the encircled spins may change state; in case I, the 
boundary spins will follow the opinion of the left domain of 
up spins which will grow. For case II, the down spin between 
the two up spins will fiip irrespectative of the size of the neig- 
bouring domains. 



boundary, is most exposed to the competition between 
opposing pressures and gives in to the larger one. This is 
what happens in case I shown in Fig[TJ The interaction 
in case II on the other hand implies that it is difficult to 
stick to one's opinion if the entire neighbourhood opposes 
it. 

Defining the dynamics in this way, one immediately no- 
tices that case II corresponds to what would happen for 
spins in a nearest neighbour ferromagnetic Ising model 
(FIM) in which the dynamics at zero temperature is 
simply an energy minimisation scheme. However, the 
boundary spin in the FIM behaves differently in case I; 
it may or may not flip as the energy remains same. In 
the present model, the dynamics is deterministic even for 
the boundary spins (barring the rare instance when the 
two neighbourhoods have the same size in which case the 
individual changes state with fifty percent probability). 

In this model, the important condition of changing 
one's opinion is the size of the neighbouring domains 
which is not fixed either in time or space. This is the 
unique feature of this model, and to the best of our knowl- 
edge such a condition has not been considered earlier. In 
the most familiar models of opinion dynamics like the 
Sznajd model Q or the voter model 4^, one takes the ef- 
fect of nearest neighbours within a given radius and even 
in the case of models defined on networks (3 , the influ- 
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encing neighbours may be nonlocal but always fixed in 
identity. 

We have done Monte Carlo simulations to study the 
dynamical evolution of the proposed model from a given 
initial state. With a system of N spins representing indi- 
viduals, at each step, one spin is selected at random and 
its state updated. After N such steps, one Monte Carlo 
time step is said to be completed. 

If iV+ is the number of people of a particular opinion 
(up spin) and N_ is the number of people of opposite 
opinion (down spin), the order parameter is defined as 
m = |A^+ — N-\/N. This is identical to the (absolute 
value of) magnetisation in the Ising model. 

Starting from a random initial configuration, the dy- 
namics in Model I leads to a final state with m — 1, 
i.e., a homogeneous state where all individuals have the 
same opinion. It is not difhcult to understand this result; 
in absence of any fluctuation, the dominating neighbour- 
hood (domain) simply grows in size ultimately spanning 
the entire system. 

We have studied the dynamical behaviour of the frac- 
tion of domain walls D and the order parameter m as the 
system evolves to the homogeneous state. We observe 
that the behaviour of D{t) and m{t) are consistent with 
the usual scaling behaviour found in coarsening phenom- 
ena; D{t) cx t-i/^ with z = 1.00 ± 0.01 and m(t) cx t^/^^ 
with z = 0.99 ± 0.01. These variations are shown in Fig. 

m 




FIG. 2; (Color online) Growth of order parameter m with 
time for two diflterent system sizes along with a straight line 
(slope 0.51) shown in a log-log plot. Inset shows the decay of 
fraction of domain wall D with time. 

We have also calculated the persistence probability 
that a person has not change his/her opinion up to time 
t. Persistence, which in general is the probability that 
a fluctuating nonequilibrium field does not change sign 
upto time t, shows a power law decay behaviour in many 
physical phenomena, i.e., P{t) oc t^^ , where 6 is the per- 
sistence exponent. In that case, one can use the finite 
size scaling relation d, 01 

P(t,L)^t-'f{L/t^/^). (1) 



For finite systems, the persistence probability saturates 
at a value oc L~" at large times. Therefore, for a; << 1 , 
f{x) cx x^" with a — z9. For large a;, f{x) is a constant. 
Thus one can obtain estimates for both z and 9 using the 
above scaling form. 

In the present model the persistence probability does 
show a power law decay with 9 = 0.235 ±0.003, while the 
finite size scaling analysis made according to H]) suggests 
a z value 1.04±0.01 (Fig. [3]). Thus we find that the values 
of z from the three different calculations are consistent 
and conclude that the dynamic exponent z = 1.02 ±0.02. 

It is important to note that both the exponents z and 
9 are novel in the sense that they are quite different from 
those of the one dimensional Ising model Q and other 
opinion / voter dynamics models [3 EOi [HI ■ Specifically 
in the Ising model, z = 2 and 9 = 0.375 and for the 
Sznajd model the persistence exponent is equal to that 
of the Ising model. This shows that the present model 
belongs to an entirely new dynamical class. 
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FIG. 3: (Color online) The collapse of scaled persistence prob- 
ability versus scaled time using 6 = 0.235 and z = 1.04 is 
shown for different system sizes. Inset shows the unsealed 
data. 

The Model I described so far has no fluctuation. Fluc- 
tuations or disorder can be introduced in several ways. 
We adopt a realistic outlook: since every individual is 
not expected to succumb to social pressure, we modify 
Model I by introducing a parameter p called rigidity co- 
efficient which denotes the probability that people are 
completely rigid and never change their opinions. The 
modified model will be called Model II in which there 
are pN rigid individuals (chosen randomly at time t = 0), 
who retain their initial state throughout the time evolu- 
tion. Thus the disorder is quenched in nature. The limit 
p = 1 corresponds to the unrealistic noninteracting case 
when no time evolution takes place; p = 1 is in fact a triv- 
ial fixed point. For other values of p, the system evolves 
to a equilibrium state. 

The time evolution changes drastically in nature with 
the introduction of p. All the dynamical variables like 
order parameter, fraction of domain wall and persistence 
attain a saturation value at a rate which increases with 
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p. Power law variation with time can only be observed 
for p < 0.01 with the exponent values same as those for 
p = 0. The saturation or equilibrium values on the other 
hand show the following behaviour: 

m, cx N-'^^p^'^^ 

Ps=a + bpf^' (2) 

where in the last equation a(~ 10~^) and /?3(~ 0.36) are 
dependent on N. The values of the exponents are ai — 
0.500 ± 0.002, /3i = 0.513 ± 0.010, /32 = 0.96±0.01. (Figs 
2] and [51 ) The variation of with p is strictly speaking 
not valid for extremely small values of p. However, at 
such small values of p, it is difficult to obtain the exact 
form of the variation numerically. 
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FIG. 5: (Color online) Scaled satuaration value of decays 
with the rigidity coefficient p. nset shows the unsealed data. 
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FIG. 4: (Color online) Satuaration values of fraction of do- 
main walls Ds and persistence probability Ps (shown in in- 
set) increase with rigidity coefficient p in a power law manner. 
There is no system size dependence for both the quantities. 

It can be naively assumed that the Np rigid individ- 
uals will dominantly appear at the domain boundaries 
such that in the first order approximation (for a fixed 
population), D oc 1/p. This would give m cx l/^/p indi- 
cating f3i = 0.5 and /32 = 1- The numerically obtained 
values are in fact quite close to these estimates. 

The results obtained for Model II can be explained in 
the following way: with p ^ 0, the domains cannot grow 
freely and domains with both kinds of opinions survive 
making the equilibrium mg less than unity. Thus the so- 
ciety becomes heterogeneous for any p > when people 
do not follow the same opinion any longer. The variation 
of rris with N shows that mg — s- in the thermodynamic 
limit for p > 0. Thus not only does the society become 
heterogenous at the onset of p, it goes to a completely 
disordered state analogous to the paramagnetic state in 
magnetic systems. Thus one may conclude that a phase 
transition from a ordered state with m = 1 to a dis- 
rodered state (m = 0) takes place for p — 0^. It may 
be recalled here that m = at the trivial fixed point 
p = 1 and therefore the system flows to the p = 1 fixed 



point for any nonzero value of p indicating that p = 1 is 
a stable fixed point. 

That the saturation values of the fraction of domain 
walls do not show system size dependence for p = 0''" 
further supports the fact that the phase transition occurs 
at p = 0. 

The effect of the parameter p is therefore very similar 
to thermal fluctuations in the Ising chain, which drives 
the latter to a disordered state for any non-zero temper- 
ature, p — 1 being comparable to infinite temperature. 
However, the role of the rigid individuals is more simi- 
lar to domain walls which are pinned rather than thermal 
fluctuations. In fact, the Ising model will have dynamical 
evolution even at very high temperatures while in Model 
II, the dynamical evolution becomes slower with p, ul- 
timately stopping altogether at p = 1. This is reflected 
in the scaling of the various thermodynamical quantities 
with p, e.g., the order parameter shows a power law scal- 
ing above the transition point. 




FIG. 6: The flow lines in the p~ p plane; Any non-zero value 
of p with p = drives the system to the disordered fixed point 
p = 1. Any nonzero value of p drives it to the ordered state 
(p = 1, which is a line of fixed points) for all values of p. 

Since the role of p is similar to domain wall pinning. 
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one can introduce a depinning probability factor /i which 
in this system represents the probabihty for rigid individ- 
uals to become non-rigid during each Monte Carlo step, 
/i relaxes the rigidity criterion in an annealed manner in 
the sense that the identity of the individuals who become 
non-rigid is not fixed (in time). If /i = 1, one gets back 
Model I (identical to Model II with p—0) whatever be the 
value of p, and therefore p — I signifies a line of (Model 
I) fixed points, where the dynamics leads the system to 
a homogeneous state. 

With the introduction of /i, one has effectively a lesser 
fraction p' of rigid people in the society, where 

p' = p{l-f,). (3) 

The difference from Model II is, of course, that this ef- 
fective fraction of rigid individuals is not fixed in identity 
(over time). Thus when p ^ 0,/x ^ 0, we have a system 
in which there are both quenched and annealed disorder. 
It is observed that for any nonzero value of p, the system 
once again evolves to a homogeneous state {m = 1) for 
all values of p. Moreover, the dynamic behaviour is also 
same as Model I with the exponent z and 9 having iden- 
tical values. This shows that the nature of randomness 
is crucial as one cannot simply replace a system with pa- 
rameters {p ^ 0, /i ^ 0} by one with only quenched ran- 
domness {p' ^ 0, p' = 0} as in the latter case one would 
end up with a heterogeneous society. We therefore con- 
clude that the annealed disorder wins over the quenched 
disorder; p effectively drives the system to the p = I fixed 
point for any value of p. This is shown schematically in 



a flow diagram (Fig. [6]). It is worth remarking that it 
looks very similar to the flow diagram of the one dimen- 
sional Ising model with nearest neighbour interactions in 
a longitudinal field and finite temperature. 

In summary, we have proposed a model of opinion dy- 
namics in which the social pressure is quantified in terms 
of the size of domains having same opinion. In the sim- 
plest form, the model has no disorder and self organises to 
a homogeous state in which the entire population has the 
same opinion. This simple model exhibits novel coarsen- 
ing exponents. With disorder, the model undergoes a 
phase transition from a homogeneous society (with order 
parameter equal to one) to a heterogeneous one which 
is fully disordered in the sense that no consensus can 
be reached as the order parameter goes to zero in the 
thermodynamic limit. With both quenched and annealed 
randomness present in the system, the annealed random- 
ness is observed to drive the system to a homogeneous 
state for any amount of the quenched randomness. 

Many open questions still remain regarding Models I 
and II, the behaviour in higher dimensions being one of 
them. In fact, full understanding of the phase transition 
occurring in Model II reported here is an important issue: 
although the phase transition has similarities with the 
one dimensional Ising model, there are some distinctive 
features which should be studied in more detail. 
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